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MODEL-BASED CONTROL SYSTEMS AND METHODS 
FOR GAS TURBINE ENGINES 


BACKGROUND OF THE INVENTION 

[0001] This invention was made with Government support under government 

contract no. F33615-98-C-2901 awarded by the U.S. Department of Defense to 
General Electric Corporation. The Government has certain rights in the invention, 
including a paid-up license and the right, in limited circumstances, to require the 
owner of any patent issuing in this invention to license others on reasonable terms. 

[0002] The present invention relates generally to systems and methods for 

controlling a gas turbine engine. More specifically, the present invention relates to 
adaptive model-based control systems and methods that maximize capability after 
deterioration, fault, failure or damage to one or more engine components or systems 
so that engine performance and/or operability can be optimized. 

[0003] Mechanical and electrical parts and/or systems can deteriorate, fail or 

be damaged. Any component in a gas turbine system, including engine components, 
sensors, actuators, or any of the engine subsystems, is susceptible to degradation, 
failure or damage that causes the engine to move away from nominal conditions. The 
effect that these upsets have on the gas turbine performance ranges from no effect 
(e.g., possibly due to a single failed sensor in a multi-sensor system) to a total loss of 
engine power or thrust control (e.g., for a failed actuator or damaged engine 
component). Control systems of gas turbine engines may be provided to detect such 
effects or the cause of such effects and attempt to compensate. 

[0004] Currently, gas turbine systems rely on sensor-based control systems, in 

which operating goals and limits are specified and controlled in terms of available 
sensed parameters. Online engine health management is typically limited to sensor 
failure detection (e.g., range and rate checks), actuator position feedback errors, and 
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some selected system anomaly checks, such as stall detection, rotor overspeed, and 
other such indications of loss of power or thrust control. When an engine component 
or system fails or deteriorates, control of the component/system is handled on an 
individual basis (i.e., each component/system is controlled by its own control 
regulator or heuristic open-loop logic), 

[0005] It is believed that presently no adequate adaptive model-based control 

systems and methods are available. 

SUMMARY OF THE INVENTION 

[0006] One embodiment of the invention relates to a method of controlling a 

gas turbine engine. The engine has sensors to detect one or more parameters and 
actuators adapted to respond to commands. The method includes receiving data from 
the sensors of the engine for one or more measured or sensed parameters, estimating a 
state of the engine by estimating one or more unmeasured or unsensed parameters 
using the data from the sensors and a predictive model of the engine, generating 
commands for the actuators based on the state using an optimization algorithm, and 
transmitting the commands to the engine. 

[0007] Another embodiment of the invention relates to a system for 

controlling a gas turbine engine, the engine having sensors to detect one or more 
parameters and actuators adapted to respond to commands. The system includes a 
state estimator adapted to estimate a state of the engine by estimating one or more 
unmeasured or unsensed parameters using data from the sensors of the engine for one 
or more measured or sensed parameters. The estimator includes a model of the 
engine. The system also includes a control module adapted to generate commands for 
the actuators based on the state. The control module includes an optimization 
algorithm for determining the commands. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0008] Figure 1 is a schematic diagram showing the layout of an engine that 

may be controlled by a system or method according to an embodiment of the 
invention; 

[0009] Figure 2 is a diagram illustrating the concept of receding horizon 

control implemented in an embodiment of the present invention; and 

[0010] Figure 3 is a schematic illustration of a control arrangement according 

to an embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0011] Embodiments of the present invention provide control systems and 

methods wherein the models, optimizations, objective functions, constraints and/or 
parameters in the control system modify, update and/or reconfigure themselves 
whenever any engine component or system moves away from nominal so that as 
much performance and/or operability as possible can be regained. Further, systems 
and methods according to embodiments of the present invention provide that the 
control system updates itself in real-time. The systems and methods may be 
automated using a computer. Embodiments of the present invention may take 
information about detected deterioration, faults, failures and damage and incorporate 
such information into the proper models, optimizations, objective functions, 
constraints and/or parameters in the control system to allow the control system to take 
optimized action given the current engine condition. Such systems and methods may 
allow any level of deterioration, faults, failures or damage to be accommodated, and 
not just deterioration, faults, failures or damage that have a priori solutions already 
progranmied into the system. Furthermore, embodiments of the present invention 
may be capable of being used to control gas turbines, such as the gas turbines in an 
aircraft engine, power plant, marine propulsion, or industrial application. 
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[0012] Figure 1 illustrates a schematic of a layout of an engine 10 as well as 

the station designations, sensors, and actuators for the engine 10. The engine 10 is an 
aerodynamically coupled, dual rotor machine wherein a low-pressure rotor system 
(fan and low-pressure turbine) is mechanically independent of a high-pressure (core 
engine) system. Air entering the inlet is compressed by the fan and then split into two 
concentric streams. One of these streams then enters the high-pressure compressor 
and proceeds through the main engine combustor, high-pressure turbine, and low- 
pressure turbine. The other stream is directed through an annular duct and then 
recombined with the core flow, downstream of the low-pressure turbine, by means of 
a convoluted chute device. The combined streams then enter the augmenter to a 
convergent-divergent, variable area exhaust nozzle where the flow is pressurized, 
expanded and accelerated rearward into the atmosphere, thereby generating thmst. 

[0013] The various actuators of the engine 10 are controlled through actuation 

inputs from a controller, such as the model predictive controller described below with 
reference to Figure 3. The various sensors provide measured or sensed values of 
parameters for monitoring and use by one or more systems. For example, the sensed 
and measured values may be used to estimate values of unsensed and unmeasured 
parameters using a state estimator, as described below with reference to Figure 3. 

[0014] It will be understood by those skilled in the art that the disclosed 

embodiments may be applicable to a variety of systems and are not limited to engines 
similar to that illustrated in Figure 1 . 

[0015] During normal operation, such engines can experience large variations 

in operating parameters, such as ambient temperature, pressure, Mach number and 
power output level. For each of these variations, the change in engine dynamics 
includes a significant nonlinear component. A control system and method for such an 
engine must adapt to such non-linear changes. 

[0016] Control systems adapted to provide control of such engines have been 

described in U.S. Patent Application Serial No. 10/306,433, GE Dkt. No. 124447, 
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entitled "METHODS AND APPARATUS FOR MODEL PREDICTIVE CONTROL 
OF AIRCRAFT GAS TURBINE ENGINES/' filed November 27, 2002, and U.S. 
Patent Application Serial No. 10/293,078, GE Dkt. No. 126067, entitled "ADAPTIVE 
MODEL-BASED CONTROL SYSTEMS AND METHODS FOR CONTROLLING 
A GAS TURBINE," filed November 13, 2002, each of which is incorporated herein 
by reference in its entirety. 

[0017] A nonlinear model predictive control (NMPC) algorithm can explicitly 

handle relevant aircraft engine control issues in a single formulation. NMPC is a 
nonlinear, multi-input, multi-output algorithm capable of handling both input and 
output constraints. Embodiments of the present invention use a dynamic model of the 
system to determine the response of the engine to control inputs over a future time 
horizon. The control actions are determined by a constrained online optimization of 
these future responses, as described in detail below with reference to Figure 3. 

[0018] The concept of receding horizon control 20 is illustrated in Figure 2. 

At time k 21, the input variables 22 {u(k), u(k'hJ), u(k-\-p-l)) are selected to 
optimize a performance criterion over the prediction horizon 23 (p). Of the computed 
optimal control moves, only the values for the first sample (u(k)), are actually 
implemented. Before the next time interval 24, 24' and calculation of another p input 
value (i.e., at u(k+l), u(k'\-2), u(k+p)), the initial state is re-estimated from output 
measurements. This causes the seemingly open-loop strategy to actually implement a 
closed-loop control. For further details, reference may be made to J.M. Maciejowski, 
Predictive Control with Constraints , Prentice-Hall London, 2002. 

[0019] Figure 3 illustrates a control arrangement implementing NMPC 

according to an embodiment of the invention. The control system 100 is adapted to 
monitor and control the physical engine plant 110 to provide substantially optimal 
performance under nominal, off-nominal and failure conditions, for example. 
"Optimal performance" may refer to different qualities under different conditions. 
For example, under normal flight, optimal performance may refer to maximizing fuel 
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efficiency, while under a failure condition, optimal performance may refer to 
maximizing operability of the engine through maximum thrust. 

[0020] The plant 110 includes sensors which sense or measure values Y of 

certain parameters. These parameters may include, for example, fan speed, pressures 
and pressure ratios, and temperatures. The plant also includes a plurality of actuators 
which are controlled by command inputs U. The plant may be similar to the engine 
illustrated in Figure 1, for example. 

[0021] The values Y of the sensed or measured parameters are provided to a 

state estimator 120. The use of NMPC requires that values of all states must be 
available. This is required since NMPC is a full state feedback controller. 
Availability of sensed or measured data is generally limited due to lack of sensors. To 
accommodate the requirements of the NMPC, embodiments of the present invention 
implement an Extended Kalman Filter (EKF) for estimating values of unmeasured or 
unsensed parameters. The EKF is described below in greater detail. 

[0022] The state estimator 120 includes a model 130 of the plant 110. The 

model 130 is used by the state estimator 120 to generate state parameters which 
include estimates of unmeasured and unsensed parameters. In a particular 
embodiment, the model 130 is a simplified real-time model (SRTM), described in 
further detail below. The SRTM is a non-linear model that can be linearized for use 
by the state estimator to determine Kalman gain values. 

[0023] The state parameters from the state estimator 120 are transmitted to a 

model-based predictive control module 140. The control module 140 uses the state 
parameters to perform an optimization to determine conmiands for the actuators of the 
plant 110. In this regard, the control module 140 includes an optimizer ISO and a 
model 160. The model 160 may be identical to the model 130 in the state estimator 
120. In a particular embodiment, both models 130, 160 are the SRTM. Using the 
SRTM, rather than a detailed, physics-based model allows the optimization to 
converge rapidly. In a particular embodiment, the optimizer 150 includes a quadratic 
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programming algorithm to optimize an objective function under given constraints. 
The optimizer determines the optimum values of control variables (i.e., actuator 
commands), and allows constraints to be specified relating to certain engine 
parameters, such as maximum temperatures, altitude and Mach number, while 
maximizing or minimizing an objective function, such as fuel efficiency or thrust. It 
is noted that, in an embodiment of the invention, constraints and objective function 
may include any of the state parameters, whether sensed, measured, unsensed or 
unmeasured. An exemplary formulation for the optimizer is described below. 


MODEL 

[0024] Physics-based, component-level models (CLM) have been employed 

for various applications, including certain control systems. A CLM is generally a 
complicated, iterative model. Such a model may require extensive processing when 
used in an optimizer, for example, thereby delaying convergence of the optimization. 
In this regard, a non-linear, analytic and non-iterating model may be implemented. 
This model is referred to herein as a simple real-time model (SRTM). In a particular 
embodiment, this model is used in both the state estimator 120 and the control module 
140. 

[0025] An exemplary SRTM for implementation of embodiments of the 

present invention has inputs including the 1) fuel flow demand, 2) exhaust nozzle area 
demand, 3) altitude, 4) Mach, and 5) delta from ambient temperature. The first two 
inputs correspond to actuator conmiands U from the control module 140, while the 
remaining three correspond to measured or sensed outputs Y from the plant 110. It is 
noted that these inputs are only exemplary and that other combinations of inputs are 
contemplated within the scope of the invention. 

[0026] The outputs of the SRTM include estimates for certain unmeasured and 

unsensed parameters. These output parameters may include core speed, fan speed. 
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fan inlet pressure, fan exit pressure, compressor inlet pressure, compressor discharge 
static pressure, compressor discharge total pressure, fan airflow, compressor airflow, 
fan inlet temperature, compressor inlet temperature, high pressure turbine exit 
temperature, fan stall margin, core stall margin, and thrust. 

[0027] An embodiment of the SRTM model depends on tables of steady state 

data to define the steady state relationships between states and inputs and on transient 
gains to represent the transient relationships. The model may be established in the 
following manner. 

[0028] First, the dynamics of the inertias are modeled. The two main states of 

the model represent the fan and core spool inertias. The first input modeled is the 
corrected fuel flow input (wfr). Then the model is changed to account for exit area 
demand as an additional input. The steady state curves may then be generated. With 
the primary states and inputs established, other outputs and other inputs are added to 
the model. With the model structure created and all of the steady state relationships 
defined, the transient *k' parameters may then be determined through system 
identification techniques. 

[0029] The embodiment of the SRTM considers the low-pressure and high- 

pressure spool speeds as two of the energy storage components, or the states of the 
model. These speeds can change state if an unbalanced torque is applied. Simply put, 
the speed increments of the engine are the integral of the surplus torques. This is 
stated mathematically as: 

-^"yEa. (Eq.l) 

where is the spool angular acceleration, J is the number of unbalanced torques, / 
dt 

is the spool inertia, and Qj is the i* torque. The origin of the torques is based on the 
concept that if the value of an input or other state is different than what the local state 
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is expecting at steady state, then it will apply an unbalanced torque to the local state. 
Using this information and Eq. 1, this idea is expressed for LP spool speed (pcn2) , 
and the HP spool speed (pcn25) as: 

pcn2 = k2* ipcn25 - gpcn25) + kwjh2 * {wf - gwfii2) , (Eq. 2) 

pcn25 = k25 * {pcn2 - gpcn2) + kwfri25 * (wf - gwjh25) , (Eq. 3) 

where pcn2 and pcn25 are the angular acceleration of the low-pressure and high- 
pressure spools, respectively, the g parameters are based on steady state relationships, 
and the k parameters are derived from transient data. Working through Eq. 2, each of 
the terms is described as follows: 

• k2 represents the aerodynamic influence of the HP spool on the LP 
spool acceleration, 

• gpcn25 is the steady state value of pcn25 based on pcn2, 

• kwfn2 is the influence of a change in wf on the LP spool acceleration, 

• gwfn2 is the steady state value of wf based on the value of pcn2. 
[0030] Similarly for Eq. 3, 

• k25 represents the aerodynamic influence of the LP spool on the HP 
spool acceleration, 

• gpcn2 is the steady state value of pcn2 based on pcn2S, 

• kwfh2S is the influence of a change in wf on the HP spool acceleration, 

• gwfn25 is the steady state value of wf based on the value of pcn25. 
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[0031] 


The two control outputs from the control module are fuel flow demand 


and exhaust nozzle area demand. The engine model inputs are fuel flow and exit area. 
Between the commands from the control and the physical inputs to the engine are the 
inner-loop control algorithm and the actuators. The models of the inner loop controls 
and actuator dynamics for both the fuel metering valve and exhaust nozzle are 
created. 

[0032] As noted above, in a particular embodiment, the SRTM is used as a 

predictive model in both the state estimator and the control module. The state 
estimator of one embodiment is an Extended ICalman Filter (EICF) that uses the 
SRTM in its nonlinear form (described above) for the time update calculation. A 
linearized version of the SRTM is used by the EKF for the Kalman gain calculation. 
Similarly, the control module of one embodiment, using a quadratic progranmiing 
algorithm, depends on a linear SRTM model to define the relationships between 
future control actions and future engine responses. 

[0033] A linearized version of the SRTM is obtained as follows. The SRTM 

can be described in general as a nonlinear ordinary differential equation (ODE): 


with the states JC/ and the inputs Ut. Taylor's theorem is used to linearize the solution 
about the current (jc, , ) value. Introducing the deviation variables (x, , ) , 


(Eq. 4) 


(Eq. 5) 


yields the following standard Taylor's expansion for the ODE in Eq. 4: 



(Eq. 6) 
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[0034] This ODE describes how the solution jc, evolves with control m, in 

comparison with the nominal solution from control . The linearized system is 
then represented by: 

M +/(x,,M;)-i. (Eq.7) 


^ a/ 

' dx 


X. + 


[0035] In the above ODE, i = 0 , since is a constant denoting the current 

value of the states. 

[0036] Moreover, for linearization about steady-state equilibrium 

solutions, /(J, m) = 0, and thus, there is no additive term f{x,,u,). However, when 

linearizing about an arbitrary current point (3c,,m,) , this additive term is a non-zero 
term that is constant over the timeframe of evolution of the linearized system. 

[0037] In addition to the ODE system in Eq. 4 that describes the dynamics of 

the system, we also linearize the output relations for both the measured outputs Zt and 
the controlled outputs yt : 


y, =K(x,,u,) 


(Eq. 8) 


using a similar Taylor's series expansion about the current values: 


y, = y,-y, 


x,+ 


du 


X. + 


M, 


(Eq.9) 
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[0038] Using the above, the identity y = h(x,,u,), and the following 

substitutions: 


A, = 


C = 


dx 


dx 


du 


(Eq. 10) 


the linear model is derived: 


(Eq. 11) 


[0039] Finally, when the control solution M,+yis determined it should be 

interpreted as additive to the constant current input, so that w,^.^ = + m,^^ . 

[0040] It is important to note the /term in Eq. 11. This term represents the 

free response of the plant. 

[0041] Thus, the embodiment of the SRTM provides simplified model that 

provides accurate and rapid convergence of the optimization. The model can be 
linearized for certain purposes. 

[0042] The linear model in Eq. 1 1 is discretized in time using the sample time 

Ts to obtain the linear discrete-time model: 


jc,^.j = Ax^ + Bu^ -f F 
y,=Cx,+Du, 


(Eq. 12) 


where 


A = / + A,r,, B = B,7;, F = f[x,^u,)T^ . (Eq. 13) 
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EXTENDED KALMAN FILTER 

[0043] The state estimator implemented in a particular embodiment is an 

Extended Kalman Filter (EKF). The EKF is a nonlinear state estimator, which is 
based on a dynamical system model. While the model underpinning the EKF is 
nonlinear, the recursion is based on a linear gain computed from the parameters of the 
linearized SRTM model. Thus the design concepts inherit much from the realm of 
Kalman Filtering. 

[0044] The EKF need not provide the truly optimal state estimate to the 

controller in order to operate adequately well. It is usually a suboptimal nonlinear 
filter in any case. However, its role in providing the state estimates to the NMPC for 
correct initialization is a key feature of the control module. 


[0045] For the EKF analysis, the SRTM is described by: 


(Eq. 14) 


where the measurement y arrives at every At seconds and the white noise variables w 
and V represent the process and measurement noises, respectively. This is a 
continuous-time dynamical system with discrete-time (sampled) measurements. 

[0046] The EKF equations can be written in predictor-corrector form. For the 
state estimation case, the predictor or time-update equations using Euler integration to 
move from continuous to discrete time are: 


* * *\ (Eq. 15) 
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where jCjt+, is a priori to the measurement step state estimate, P^^, is the a priori 
estimate error covariance, W is the discrete-time process noise covariance (after 
scaling by At), and Af^ is the discrete-time transition of the linearized system, or: 


A, =/ + A,r, (Eq. 16) 

[0047] The linear discrete time measurement matrix C is defined as: 


C»=^(i;...P*."*)- (Eq.l7) 


[0048] Next the Kalman filter gain is computed using: 

K = P;,,C,[r + C,P;,,C,Y (Eq. 18) 


[0049] The corrector or measurement update equations are: 


-k(r^-c,p:,,c,Yk\ 


P = P" 


(Eq. 19) 


OPTIMIZER FORMULATION 


[0050] Embodiments of the control module include an optimizer adapted to 
maximize or minimize an objective function while satisfying a given set of 
constraints. In one embodiment, the optimizer uses a quadratic programming 
algorithm. As described above, the control module uses a dynamic model of the plant 
to perform simulations over a specific horizon, the model being the SRTM in one 
embodiment. 
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[0051] In an exemplary embodiment, the control module is designed to control the 
fan speed PCN2R, and the pressure ratio DPP (yif=PCN2R, y2t=DPP) using the 
combustor fuel flov^/ , finvdmd, and the afterburner, aSxdmi as the manipulated inputs 
{uit-fmvdmd, U2t = aSxdmi), subject to magnitude and slew rate constraints imposed 
by the hardware limits on the actuators for the two manipulated inputs. In addition to 
these constraints, the optimization in the control module will also be performed 
subject to other operational/safety constraints like stall margin, combustor blowout, 
maximum T4B, minimum and maximum PS3, maximum N25. 

[0052] A quadratic programming (QP) algorithm with the linearized dynamic 
model is used along with a quadratic objective function and linear constraints. The 
QP problem is convex and can be solved readily with available QP software. 
Moreover, since the linearization is performed repeatedly at each time sample about 
the corresponding operating point, it accounts for the nonlinearities encountered 
during dynamic transients over the flight envelope. To implement the control module 
using the QP formulation, the nonlinear SRTM is linearized about the current state 
estimate obtained by the EKF and the current inputs w,_i , and then discretized in 
time using the sample time T^, to obtain the linear discrete-time model. 

[0053] The main control objective is to track changes in the references for the two 
controlled outputs. The optimization objective function is postulated as a standard 
quadratic function to be minimized over a future prediction horizon n/i, using the 
piecewise constant inputs over a future control horizon nc. More specifically, the 
objective function to be minimized is: 


1 > 


Q'^y..^ 


t+nh 


1 \ 


(Eq.20) 
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where Ay,+, = y^j+i - denotes the error between the output reference and the 
predicted output at a future sample r+i, Am,^.. = u^^. - m,^.._i = Am,^.,. = u^^. - m^^._, 

denotes the change in the manipulated inputs at sample t-hi relative to the value of the 
inputs at the previous sample, and Q and R are symmetric positive definite weighting 
matrices. The weighting matrices Q and R and the prediction and control horizons nh, 
nc, respectively are tuned for optimal performance and stability. The objective 
function is to be minimized as a function of the future control action values 
= l,...,nc assuming that = ,/ = nc + l,...,n/i . The objective function 
is calculated over the future prediction horizon nh using a linear discrete time model. 

[0054] The predicted values of the outputs y^+i in terms of deviations from the 
current measured value y, are given by the following relation: 


, where =(/„,®C) 





-^1+1 






_yt+nh _ 




A 


X. + 


B 0 •• 0 
AB B 
: : ■•. 0 


U 


t+nh_ 


1=1 

nh 


(Eq. 21) 


0 1 


[0055] Note that by its definition, 3c, = 0 , which is utilized in the above relation to 
obtain the predicted outputs over the prediction horizon nhj as a linear function of 
the future control action over the control horizon nc. Moreover, the changes in the 
control inputs Am,^. = m,^. -m,+,_i are denoted by the following compact relation: 
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Am, 


" 1 

0 

... 0" 




-1 

1 

0 0 




0 


•. 0 




0 

0 

-1 1 



(Eq.22) 


[0056] Using the above relations, the predicted value of the objective function to 
be minimized over the prediction horizon is given by the compact relation: 


(Eq.23) 


[0057] The above quadratic objective function is to be minimized with respect to 
the future control moves , subject to all input and output constraints. In particular, 
the input constraints consist of the min/max magnitude and rate of change constraints: 


"min ^^t-l ^«max 


(Eq.24) 


[0058] In addition to the above input constraints, which are typically hard 
constraints, there may be other state/output operational/safety constraints (e.g., 
minimum stall margin, maximum core speed, combustor blowout). In one 
formulation of the NMPC, a logic to generate the output reference trajectory and 
update the constraints for changes in the control actions (fuel flow and A8) is used to 
enforce these operational/safety constraints. However, it is possible to enforce these 
operational/safety constraints directly using a linear model for the prediction of the 
relevant state/output variables over the prediction horizon. For instance, in order to 
enforce the maximum limit on the core speed, which is a measured variable and the 
2°^ state in the SRTM, the constraint can be accounted for using the linear discrete- 
time model: 


•^2,/+/ — -^2 max ^24 ^2 


(Eq.25) 
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[0059] Note that, unlike the input constraints, these state/output constraints rely 
on the model predictions and thus are subject to plant-model mismatch over the 
prediction horizon. Thus, to avoid potential infeasibility, these constraints are 
typically included as a soft constraint. Thus, the overall QP problem to be solved at 
each time sample for the NMPC is given below: 


min J^^ = UlH.U, + HJU, +Wfi 
subject to the constraints : 

U <U <U (Eq. 26) 


[0060] In the above QP formulation, the constant term Ho in the quadratic 
objective function is ignored, denotes the violation in the soft, output/state 

constraints, W is the penalty on the soft constraints, YJ^ and YJ^^ are the minimum 
and maximum limits on these output/state constraints in terms of deviations from the 
current values and L^U^ + denotes the predicted values of these output/state 
constraints over the prediction horizon using the linear discrete time model. 

[0061] The solution of the QP problem in Eq. 26 yields the optimal control 

trajectory U ^ over the control horizon nc. The optimal values for the first sample, i.e. 

, yields the absolute value of the control action, = m^_i + . This optimal control 

input is implemented and the QP problem is updated and solved at the next sample 
along with the EKF. 

[0062] The quadratic-progranrmiing based optimizer and control module rely on 
the predictions of the engine variables over the future prediction horizon. In the 
presence of a plant-model mismatch, the model predictions used in the control module 
can be incorrect and can lead to controller performance degradation or even 
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instability. In embodiments of the present invention, the plant-model mismatch is 
addressed by including a corrective term on the model used for the prediction. In 


provides the mismatch between the current output measurements y, and the model 
predictions for these outputs y, -h{x~^^,p^,u^) with the current state estimates jc, . 

The current value of this feedback correction term can be used as a constant 
correction term in the linearized discrete-time model. More specifically, this constant 
term can be included in the constant vector F to obtain the corrected linear model that 
can be used for prediction: 


[0063] The above correction term that accounts for the mismatch between the 
measured and predicted outputs, along with the fact that the quadratic objective 
function formulation in terms of deviations in the control actions effectively amounts 
to an integral action with respect to the error between the output reference and the 
predicted outputs, allows an offset-less control even in the presence of plant-model 
mismatch. 

[0064] In another formulation of the QP problem, an infinite prediction horizon 
may be implemented. In this regard, the control objective function is extended to an 
infinite prediction horizon, resulting in positive impacts on stability and robustness. 
The penalty for using the infinite prediction horizon is an increase in the 
computational cost to achieve a solution. To counter this penalty, a compact and 
efficient calculation of the infinite horizon term has been developed. 

[0065] The standard quadratic objective function of Eq. 20 involves a quadratic 
cost on the tracking error Ay^^, over a prediction horizon n^, and a quadratic cost on 

the control action over a control horizon (n^ « ), where it is assumed that 

the control action is constant after the control horizon, i.e. = u^^^ = "i+m:+i = ' 


particular, at each time sample t, the term 



(Eq. 27) 
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A larger control horizon enables improved control performance, however the 
optimization problem and hence the computational burden grows with the control 
horizon, thereby limiting the control horizon due to real-time implementation issues. 
On the other hand, a larger prediction horizon enables improved stability and 
robustness, hence the prediction horizon is typically chosen to be significantly larger 
than the control horizon. 

[0066] In the case of a large prediction horizon n^^ , the objective function in Eq. 

nh 

20 involves an expensive calculation of the tracking error terms ( ^AyJ^iQ^t^i) 

l=/IC+l 

beyond the control horizon . This increases the computational burden and limits 

the choice of the prediction horizon due to real-time implementation issues. The 

use of infinite prediction horizon improves the stability and performance of the 
controller without adding undue computational burden. A significantly more efficient 
alternative is proposed to evaluating the quadratic cost due to the tracking error over 
an "infinite" prediction horizon with minimal computational overhead. 

[0067] In particular, consider the quadratic objective function over an infinite 
prediction horizon: 

1 r oo i=nc 1 

\ [ nc i=nc 1 1 ~ 

=^\T^yLQ^y.H + Z^L-xR^..... +t T^yI.,Q^y... (Eq. 28) 

Jnc +Jnc.-- 

[0068] Note that due to the assumption of constant control action beyond the 
control horizon (e.g., w,^.^_, =m,+^ ="r+nc+i = " the quadratic cost of the control 
action based on Am,+^+,._, is zero and omitted from the objective function. The 

objective function is factored into two terms, where the first term is the standard 
objective function Jnc corresponding to a prediction horizon nh same as the control 
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horizon nc. It is given as a quadratic function of the control action 

J^^^hl-'H^U + flU (Eq.29) 

[0069] The second term =— ^Ay,'^,.QAy,+, is the remaining quadratic cost 

on the tracking error beyond the control horizon, and needs to be computed as a 
function of the control action in a compact and efficient manner. We will henceforth 
focus on calculating this tracking error term over the infinite horizon. In fact, we will 
calculate a slightly modified term: 

i f^aAyLQ^y,,,, (Eq.30) 

with an exponentially decaying weighting factor a. given by 
^nc+i - 1» ^i+i ~ < 1) • The use of such an exponentially decaying weighting 

factor is motivated by several factors: (i) Due to modeling errors, model predictions 
over future get less accurate with increasing horizon, hence the decaying weighting 
factor reduces the weights on tracking error with increasing samples in future and 
gives more weight to tracking error in the immediate future, (ii) In some cases, one or 
more limiting constraints become active and inhibit an offset-less tracking, i.e. the 
tracking error term Ay,^.,.does not decay to zero over the infinite horizon. In such a 

case, the exponentially decaying weighting factors, (with a<l)) is necessary to 

ensure that the sum of tracking error terms over an infinite horizon is still bounded 
and can be minimized. 

[0070] The tracking error terms Ay,+, in Eq. 20 correspond to the outputs of the 
system: 
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3c, . ; . , = Ax, . , + Bu,.„, , + F 
m+i ,+nc-i ^ (Eq. 31) 


Starting from the initial state 3c,^„p and constant inputs «,+„<.-! • It is assumed that the 
above dynamic system is stable (i.e., all eigen values of A are within the unit circle, 
else the states and hence the outputs would go unbounded over the infinite prediction 
horizon). For such a stable system, the final steady state corresponding to the 
constant input is given by: 


(Eq. 32) 


or. 


X, =(/-A)-'[fiM,,_,+F] 

y, = [C(/ - A)-' B + + C(/ - A)-' F = + ' 


(Eq. 33) 


[0071] Defining the deviation variables ic,+, = -^i. and y,+,. = yj+i -ys' 
system dynamics are given by the simplified set of equations: 


(Eq. 34) 


[0072] Thus, the infinite horizon tracking error term is given by: 
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1 *** 

CO oo 

(Eq. 35) 

[0073] Note that in the above equation, the first term is a constant, which is 
independent of the control action and can be omitted from the optimization objective. 
Moreover, the term iy^^^^i ~^u"/+/ic-i) denotes the steady state error between 
the output * references and the controlled outputs, which will be assumed to be zero. 

Also, the sunmiation term ^{af^x^^^Y QCial^x^^^) is evaluated in a compact 

l=/IC+l 

closed-form as = ^,+;,c+i^/+nc+i2^,+„c+i » where Q is a synmietric positive definite 
matrix that is the solution of the Lyapunov equation: 

Q - A^a° 'Ga° ' A = C^QC . (Eq. 36) 

[0074] Finally, x,^^^, = - 2, = - (/ - AY' (Bu^^^_, + F) , where the 

state jc,^.„^ = G^^^U + V,^^ is a function of the control inputs U and the free response 
corresponding to F, Thus, 

x„„c.x = ^..nc + ^ = AG„^U + A V„„, +F= G„^,,U + V,,^,, (Eq. 37) 

and 

= {G„^.x - [0 • • • 0 (/ - A)-' B] )U + { - (/ - A)-' F} = G,,„,„f/ + V;,„„ 

(Eq. 38) 
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[0075] Substituting these relations in Eq. 29, the following compact relation is 
obtained for the infinite horizon tracking error term: 


(Eq. 39) 


which is another quadratic expression in the control action U similar to the objective 
function One over the control horizon in Eq. 29. Thus, evaluating the matrices: 


0 



(Eq. 40) 


the infinite horizon tracking error term in Eq. 39 can be obtained in a compact and 
efficient manner. Finally, note that cc^^^ = 1 and the forgetting factor a < 1 can be 

tuned to shorten or lengthen the extent of the infinite horizon tracking error that 
contributes to the overall objective function and the rate of decay of their relative 
weighting. A larger value of a will lengthen the effective terms in the infinite horizon 
thereby increasing the stability characteristics. However, in the presence of modeling 
errors, the undue weight on distant future tracking errors will degrade the transient 
performance. A judicious tuning of the factor a will enable increased stability as well 
zs improved performance. 


[0076] The present application is related to U.S. Patent Application Serial No. 

10/306,433, GE Dkt. No. 124447, entitled "METHODS AND APPARATUS FOR 
MODEL PREDICTIVE CONTROL OF AIRCRAFT GAS TURBINE ENGINES," 
filed November 27, 2002, and U.S. Patent Application Serial No. 10/293,078, GE 
Dkt. No. 126067, entitled "ADAPTIVE MODEL-BASED CONTROL SYSTEMS 
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AND METHODS FOR CONTROLLING A GAS TURBINE," filed November 13, 
2002, each of which is incorporated herein by reference in its entirety. 

[0077] Exemplary embodiments of control systems and methods are described 

above in detail. The systems are not limited to the specific embodiments described 
herein, but rather, components of each system may be utilized independently and 
separately from other components described herein. Each system component can also 
be used in combination with other system components. 

[0078] While the invention has been described in terms of various specific 

embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the claims. 
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